Collective flow in central Au-Au collisions at 150, 250 and 400 A MeV 
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Radial collective flow and thermalization are studied in gold on gold collisions at 150, 250 and 400 
A MeV bombarding energies with a relativistically covariant formulation of a QMD code. We find 
that radial flow and "thermal" energies calculated for all the charged fragments agree reasonably 
with the experimental values. The experimental hardware filter at small angles used in the FOPI 
experiments at higher energies selects mainly the thermalized particles. 
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In recent years central-collisions studies became a focus of attention in the intermediate energy domain (100-500 
A MeV) One of the measurables concentrated on by the experiments is connected to the flow. It is a well- 

known fact [ 3| , that at larger impact parameters there is a sideward flow and a squeeze-out flow; these quantities 
were measured and calculated || . The sideward and squeeze-out flows mostly disappear in very central collisions. 
However, calculations J a] predicted a large collective energy (radial flow) in central collisions and this was confirmed 
later by experiments (€ \. This collective energy can be visualized with a blast model [00] , where the system expands 
spherically around the center of mass. 
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FIG. 1. Number of particles that collided less than three times (solid line), and without collisions (dotted line) for three 
incident energies, at b — 0.5/m impact parameter, as the function of time. Initially (t — 0) the nuclei are at a distance of 2 fm 
from each other and the total overlap for free evolution would occur at 29, 23 and 19 fm/c for the three energies, respectively. 

Recently we deduced a momentum dependent, relativistically invariant two-body force ||, which can be applied 
to QMD calculations ||. To check the validity of the force, we made detailed calculations for central gold on gold 
collisions in the energy domain 150-400 A MeV || and compared the results with experimental data EJ]. The 
agreement turned out to be highly satisfactory; even the number of the intermediate mass fragments (IMF) is very 
close to the experimentally measured value. It is worth examining what can we learn about the radial flow using this 
model. 

First we study to what extent one may speak about thermalization of the nucleons in central and semi-central 
collisions. The system is usually assumed to be thermalized if the nucleons collided a few (~ 3) times. From Fig. |l| 
one can see that for higher energies (250 — AOOAMeV) only 20% of the nucleons did not collide more than twice (and 
may be considered as unthermalized) . At low energy (150 A MeV) this fraction amounts to 60%, due to the large 
Pauli blocking (for 150 AMeV 65% of the possible collisions were blocked, while for 400 AMeV the blocking was only 
25%). We followed the path of some nucleons in collisions with 400 A MeV energy. Those nucleons which did not 
collide at all (less than 3%) are generally positioned in the outer layer of the colliding nuclei. 
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FIG. 2. Contour plots averaged over 360 events for 400 A MeV collisions with b=0.5 fm impact parameter. Shown are 
the raw proton distribution (far left), the distribution with a 9i a b > 1-2° experimental hardware filter (center left), with an 
additional 8i a b ~ 21° filter (center right), and with an extra 9i a b < 30° filter (far right) in the scaled momentum - rapidity 
plane. The hardware filter around 0; a j, ~ 5° is not visible. The contours are separated by factors 1.5. 

Assuming that nucleons having collided at least three times are thermalized, one may conclude from Fig. [l], that 
at least for high energies the thermalization rate is ~ 80%. For a spatially homogeneous system this would mean an 
isotropic momentum distribution. For this reason it is interesting to examine the momentum distribution by making 
a contour plot of the invariant cross section d 2 o~ /pTdpTaly of the outgoing single protons in the two dimensional space 
of the transverse momentum scaled with the projectile momentum, plj,, and the rapidity y scaled with the projectile 
rapidity y p in the center of mass system. We note that in the energy domain of interest (up to 600MeV) the scaled 
rapidity and scaled momentum are equal to each other within 5%. As a consequence, a thermalized distribution 
should appear as circles in the y — p plane. First we examined the contour plots without experimental filters at 
b = 0.5 fm and 400 A MeV collision energy after 150 fm/c evolution. One can see from Fig. that in contrast to 
a thermalized system, the raw data do not show an isotropic distribution (see the most inner - highest multiplicity 
contour, having a distortion factor ~ 1.5). However, applying only the seemingly negligible small laboratory forward 
angle filter (excluding particles at 6i a b < 1.2°), which means the exclusion of a small number of protons only, a 
nearly isotropic distribution is recovered. The interpretation of this result is the following: the unthermalized (not 
yet collided) protons leave the collision region with high energy in the forward direction; thus, the remaining proton 
distribution is already thermalized. Applying further angular filters does not change this behavior. 
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FIG. 3. Contour plots for proton distributions at 150 A MeV, b=1.5 fm (left block of 2 panels) and b=7.5 fm (right block of 
2 panels) in the scaled momentum - rapidity plain. The left figures of each block show the raw data, while the righthand ones 
show the data with the most relevant 8i a t > 1-2° and 9i a b < 30° hardware filters. 

Figs. H and || show contour plots for protons at 150 and 400 A MeV, respectively, for central and peripheral events. 
In order to compare our result with the experimental one, where the effect of the filters in the range 1.2° < Oi a b < 30° 
was averaged out, we took only the relevant 1.2° < 9i a b and 6i a b < 30° filters. One can see that for central events the 
filtered distributions are almost isotropic (except for the &i a b < 30° filters), however, the raw, unfiltered values are 
not. This fact is very pronounced in the case of 150 A MeV. In the case of large impact parameters the distributions 
are always distorted. The same observation was made in connection with the experimental measurements [|lj. 

Experimentally the collective radial flow is determined from the kinetic energy distribution of the large fragments. 
However, since at 400 MeV the total mass of all large fragments for central collisions is less than 3% of all the 
particles, such a method produces very poor statistics in our case. In order to determine the radial flow we considered 
to possibilities: calculating the flow in the 1. gas scenario, when all the nucleons are considered to be individual 
objects contributing to the radial flow; or in the 2. cluster scenario, when all the charged fragments (single protons 
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and clusters) are taken into account. For comparison we present results for both. 
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FIG. 4. Same as Fig. | at 400 AMeV. 

In order to determine the radial flow of the nucleons we divided the solid angle into 32 equal pieces and calculated 
the flow velocity (3^ within each sector k (fc=l,32) as 



i(fc) 



E 



/3(fe) = 
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(1) 



where u^) is the four velocity associated to sector k and the summation is extended for all the charged particles 
within the solid angle £lk, with m, being the mass of a fragment, while is the position of the charged particle at the 
end of the calculation (150 fm/c). We repeated the same for the gas scenario with each being the nucleon mass. 

The fluctuation of the flow velocity is less than 2% in all the 32 sections for each energy at b = 0.5 fm both for the 
cluster and gas algorithm. Furthermore, the fluctuation of the number of particles within the sectors is at most 3-4%. 
The extracted values are also stable against changing the end-time of the calculation: from the freeze-out time up 
to twice the freeze-out time the change of the flow energy is 10% and does not change further for larger end-times. 
These results show that the radial flow can be determined in a very reliable way. Encouraged by the isotropy of the 
system and by the lack of other collective excitations (side flow or squeeze-out) we define a disoriented ("thermal") 
energy as the rest kinetic energy[|, 
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where Ai is the mass number of a cluster, M is the total multiplicity of the clusters (and M = A, the total number 
of particles, for the gas algorithm). 

In Table | we give the calculated flow velocity, the flow energies, and the "disoriented/thermal" energies evaluated 
both for the gas and for the fragments (cluster). As a comparison, we give the values extracted from Ref. [Q, where 
they use a blast wave model fit to the experimentally measured kinetic energies of heavy fragments. We note that 
these quantities are model dependent, and evaluating them we used the natural way for the QMD model, not the 
experimental procedure. We find the agreement surprising. 

Finally we make a remark comparing our result to the one of the EOS Q group. In this experiment the selection 
criteria were based on the multiplicity only and the deduced values of the radial flow velocity are considerably 



1 The energy defined here is not exactly the thermal energy, since a thermal energy can be defined for completely thermalized 
systems only in the local rest frame. As the complicated experimental procedure obtaining the thermal energy cannot be 
reproduced within the given model, we consider the disoriented energy Edo to be close to the thermal one. 



3 



lower compared to the flow velocities of the FOPI experiment. In our model we find that the multiplicity trigger is 
not sufficient enough for selecting the central collisions; we got a considerable amount of events contributing from 
b = 4.5 — 5/m. As a result, using the multiplicity trigger, our flow velocities, (3, are reduced by ~ 15 — 20% in 
agreement with the EOS result. Consequently the "thermal" energy is increasing by the same amount. 

In conclusion, we have investigated collective radial flow in the case of Au+Au collision at 150, 250 and 400 A 
MeV. The results were compared to FOPI experiments jij and we found a reasonable agreement. Furthermore the 
experimental setup automatically filters out the unthermalized particles for higher energies. These results suggest 
that the compressibility and the momentum dependence of the used force is highly satisfactory. The application of 
our force for higher energies when particle creation is important, is in progress. 

TABLE I. Average flow energies (MeV), "disoriented/thermal" energy and flow velocity at 400 (upper), 250 (middle) and 
150 A MeV (lower part) at b=l fm impact parameter for the gas (g) and cluster (c) assumption. 
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g (400) 
c (400) 


63.2±2.5 
59.7±2.5 


56.8±6.3 


37.2±3 
32.7±2 


32.8±6.3 


0.345±0.006 
0.338±0.006 


0.334±0.017 


g (250) 
c (250) 


41.2±1.2 
34.8±1.0 


34.0±3.9 


18.2±0.5 
16.8±0.6 


21.5±3.9 


0.285±0.002 
0.264±0.002 


0.263±0.014 


g (150) 
c (150) 


21.6±0.5 
18.3±0.3 


19.9±2.3 


12.9±0.8 
10.2±0.5 


12.6±2.3 


0.21±0.002 
0.194±0.002 


0.204±0.011 
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